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ABSTRACT 

0^ \ We present spectroscopic confirmation of the "Pisces Overdensity" , also known as "Structure J", 

. a photometric overdensity of RR Lyrae stars discovered by the Sloan Digital Sky Survey (SDSS) at 

' an estimated photometric distance of ~ 85kpc. We measure radial velocities for 8 RR Lyrae stars 

\ within Pisces. We find that 5 of the 8 stars have heliocentric radial velocities within a narrow range 

{3J). of — 87kms~^ < Vr < — 67kms~^, suggesting that the photometric overdensity is mainly due to a 

^ ' physically associated system, probably a dwarf galaxy or a disrupted galaxy. Two of the remaining 3 

\ stars differ from one another by only 9 km/s, but it would be premature to identify them as a second 

system. 

Subject headings: Galaxy: halo, Galaxy: structure, galaxies: dwarf, stars: RR Lyrae 
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1. INTRODUCTION 

In the currently favored picture of galactic struc- 
ture formation, the Milky Way had a tumultuous early 
history. Continuously bombarded from an early age 
by other galaxies big and small, the Milky Way has 
been roiled by mergers — a picture first put forth by 
ISearle fc ZimJ (|1978| l. The fossil record of the Milky 
Way's history can be seen today in the debris left over 
from this cosmic carnage - the stellar streams that ap- 
pear to be the shredded remains of galaxies past. The 
SDSS has had a major impact in breaking open this field. 
A variety of photometric techniques have been devel- 
oped to search through the large photometric database 
of the SDSS to uncover these systems (e.g. Newberg et 
al. 2002, Willman et al. 2005, Belokurov et al. 2006, 
Grillmair et al. 2006). The primary search algorithms 
rely on a combination of correlations in position on the 
sky and position in the color-magnitude diagram (CMD). 
To date, a significant number of such structures have 
been located. However, the census of structures is still 
very far from complete, particularly at large Galacto- 
centric distances. The key issue is sensitivity: at very 
low intrinsic luminosities these objects can only be dis- 
entangled from foreground stars interior to '--^SOkpc and 
SDSS star counts are only sensitive to surface bright- 
nesses above ^ 3 mag/arcsec~^. A novel t echnique, 
first recognized bv lKinman fc WirtanenI (jl963l ). to push 
the limits of finding substructure within the Milky Way 
is to make use of the photometric variability of RR Lyrae 
stars. 

The repeated scans of SDSS Stripe-82 have allowed 
the discovery of Milky Way substructures via their pop- 
ulations of RR Lyrae stars, which are standard candles, 
can be seen to large distances, and show distinctive light 
curves. The structure discussed in this Letter was first 

^ This paper includes data gathered with the 6.5 meter Magellan 
Telescopes located at Las Campanas Observatory, Chile. 

^ Observatories of the Carnegie Institution of Washington, 813 
Santa Barbara Street, Pasadena, CA 91101 

^ Department of Astronomy, The Ohio State University, 4051 
McPherson Laboratory, Columbus, OH, 43210 

Department of Astromony, University of Washington, Box 
351580, Seattle, WA 98195-1580 



identified as an overdensity of RR Lyrae stars by Sesar 
et al. (2007) as part of a larger study of all RR Lyrae 
stars in Stripe-82. They termed this overdensity "Struc- 
ture J", one among several such structures located. In 
a subsequent work, Watkins et al. (2009) independently 
found what appears to be the same structure, which they 
termed the "Pisces Overdensity" . The photometric iden- 
tification of these overdensities is a crucial first step in 
locating more streams and galaxies interior to the virial 
radius of the Milky Way. However, photometry provides 
two phase-space dimensions very accurately (angular po- 
sition) and a third only very crudely (distance) . In order 
to show definitively whether photometric overdensities 
are, in fact, truly part of a common structure as opposed 
to a chance concentration, it is necessary to obtain data 
for additional phase-space dimension(s), which could in 
princi ple be either radi al velocities or proper motions 
(e.g.. ISimon et al.]l2007f ). In this Letter, we report on 
the first results of a campaign to confirm distant struc- 
tures in the Milky Way halo as defined by RR Lyrae and 
giant stars. 

2. OBSERVATIONS 

2.1. Target Selection 

In Figure[T]we show all of the identified RR Lyrae stars 
in the apparent magnitude range 19.9 < Vq < 20.8 that 
he in SDSS Stripe-82, which covers a 3° x 118° stripe 
along the celestial equator. The positions and Vq magni- 
tudes of these objects co me from the latest Stripe-82 RR 
Lyrae catalog (jSesar et a l. 2009). There is a very clear 
overdensity in the RA interval 335-360. With only about 
30% of the stripe area, this sub-region contains 24 of the 
31 RR Lyrae stars (77%). For our first campaign we fo- 
cused on the concentration located at RA~ 355° and se- 
lected 8 objects near this overdensity based on their close 
2D proximity and their similar median magnitudes. 

2.2. Light Curve Analysis 

Obtaining a single-epoch spectrum for an RR Lyrae 
star yields only a velocity of the stellar photosphere, 
which for these pulsating variables can deviate from the 
systemic velocity by several 10s of kms~^. However, 
the dense temporal sampling of Stripe-82 enables us to 
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obtain accurate ephemerides by phase folding the ~ 70 
epochs in the SDSS photometric archive. We can then 
correct velocity measurements made at known phase to 
the barycentric velocities for these stars (Joy 1938). We 
extracted the lightcurves from the SDSS archive and de- 
termined their periods based on a v ariant of the "Phas e 
Dispersion Minimization" technique (jStellingwerf lllQTSi ) . 
We first identified potential bad data points by perform- 
ing a regression on g vs. r flux, and flagging 2.5 a outliers 
(as determined from the scatter ~ not the formal errors). 
While the colors of RR Lyrae stars change during their 
pulsation cycle, the 2.5 a criterion allows an adequate 
range for normal color variation. We then removed near- 
achromatic points that were substantially fainter than 
the remaining points, which are probably due to some 
joint photometric anomaly, but in any case would not 
fit any RR Lyrae-like light curve. Finally, we varied 
the period and minimized the sum of the squares of the 
photometric differences between successive points. Our 
derived g-band lightcurves for our 8 target objects are 
shown in Figure [21 In all but one case, periods derived 
from the g-band data were nearly identical to periods de- 
rived from the r-band data. Based on these differences, 
we estimate the period error to be (j{P) = 3 x 10~^ day. 
Interestingly, our derived periods are surprisingly simi- 
lar, with a dispersion of only 0.026 days. The objects are 
all RRab type variables and with median brightnesses of 
about g ~ 20.5 these stars are at a distance of about 
85 kpc. We have also estimated metallicities using the 
periods and amplitudes for the 8 stars using equation (6) 
of Sandagc (2004j), assuming that the V-band amplitude 
is equal to the (/-band amplitude. The metallicities range 
from [Fe/H]'--^ -1.3 to -2.1. While these estimates are nec- 
essarily crude, we can conclude that all of the 8 stars are 
relatively metal-poor. 

2.3. Spectroscopic Observations 

The spectroscopic observations were obtained using the 
Magellan Echellette Spectrograph (MagE^; Marshall et 
al. 2008) mounted on the 6.5-m Clay Telescope at Las 
Campanas Observatory. In all cases two equal-length ex- 
posures bracketed an observation of a ThAr lamp. The 
data were reduced using a pipeline written by D. Kel- 
son following Kelson (2003). Post-extraction processing 
of the spectra was done with the IRAF*^ ECHELLE 
package. A 1-arcsec slit was used, resulting in a reso- 
lution of ^4100, and the signal-to-noise ratio of the final 
spectra ranged from 7 to 16 per 0.36A pixel at 4700A. 
Exposure times varied from 3600s to 4000s depending 
on the observing conditions. Velocities were measured 
with the IRAF FXCOR routine using a MagE observa- 
tion of the bl ue metal-poor star CS22 874-009 {Vheiio = 
-36.6 kms~^. fPreston fc SnedenI |2000[ ) as the template. 
The cross correlations were made on the wavelength in- 
terval 4000A - 5600A with the hydrogen lines masked 
out. 

We adopt the time-averaged velocity of the pulsation 
curve as the center-of-mass velocity of the star. Inte- 
gration of detailed velocity curves of the RRa variables 

^ http: / / www.lco.cl/telescopes-information / magellan/instruments / 
® IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under a cooperative agreement 
with the NSF. 



WY Ant, XZ Aps, DT Hya, and RV Oct (Preston 2009, 
in preparation) shows that the pulsation velocity is equal 
to the star's time-average velocity at phase 0.37, reck- 
oned relative to maximum light. Our observations were 
all made as close to this phase as possible, and velocity 
corrections were applied adopting k — 92.7 kms~^ per 
unit phase for the mean slope of the pulsation velocity 
curve at phase 0.37 for these four stars. 

A journal of the observations is presented in Table 1. 
Column one lists the heliocentric Julian date at mid- 
observation followed by the RA and Dec of the star, the 
median g-band magnitude, the adopted pulsation period, 
the observed heliocentric velocity and the error as esti- 
mated from the FXCOR measurement, the number of 
photometric datapoints in the lightcurve, the heliocen- 
tric Julian date at maximum light (zero phase), the pul- 
sation phase at mid-observation, the velocity correction, 
and the final adopted heliocentric velocity. 

3. RESULTS 

We were able to obtain accurate radial velocity mea- 
surements for 8 RR Lyrae stars associated within the 
apparent Pisces Overdensity. We show in Figure [3] a his- 
togram of the heliocentric radial velocities measured for 
our targets. There is a clear velocity peak, containing 
5 of the 8 objects, at approximately — 75kms~^. The 
other 3 stars also have an interesting velocity distribu- 
tion to which we return later. Figure [4] shows our tar- 
gets coded by velocity and magnitude. We also show 
in this figure othe r RR Lyrae stars identified in the field 
(jSesar et al.11200 9') fainter than Vq = 20. The spatial con- 
centration at (RA,DEC)- (355°, -0.3°) is suggestive of 
a coherent structure. However, our radial velocity mea- 
surements reveal a more nuanced picture. One of these 
objects (cyan square) has an apparent magnitude that 
puts it farther out in the halo than the rest of the con- 
centration. The three central stars (marked as triangles 
in the figure) within this concentration have measured 
radial velocities^ —198 km s""'^ < Vr < — 155kms~^ . 
Four stars in the central concentration and a fifth fur- 
ther away, however, have radial velocities that lie within 
the main velocity peak. We therefore seem to have two 
co-spatial velocity structures, one robustly identified in 
phase-space and one that is suggestive but too sparsely 
sampled at present to warrant detailed analysis. While 
few data exist at these distances in the stellar halo, a 
Besangon model (jRobin et al.l l2003l ) of the heliocentric 
velocity distribution at this position in the sky suggests 
that the expected velocity distribution of the smooth 
stellar halo is centered at ^ — 120kms^^ with a disper- 
sion of ~ 90kms~^. In Section [4l we focus exclusively 
on the main velocity peak. 

4. DISCUSSION 

4.1. Is the Velocity Group Physically Associated? 

In order to determine the random chance that a group- 
ing of 5 out of 8 stars could have velocities within 
20kms"' of one another, we performed a Monte Carlo 
i|§^t by random samplings from a Gaussian distribution 

^ The single object with differing g-band and r-band period 
has a barycentric velocity of — 197.7kms~^ if the g-band period 
is adopted and — 182.7kms~^ if the r-band period is adopted. 
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with dispersion of 90 km s^ . We find this probability to 
be less than 0.6%. This is a fair test (and not an a pos- 
teriori justification of a curious velocity structure found 
serendipitously) because this corresponds to exactly the 
velocity signature we were looking for when we under- 
took the observations. The group is therefore physically 
associated at high confidence in a fully Bayesian sense. 

4.2. Bound or Unbound? 

When looking at the spatial and kinematic informa- 
tion, the first question one must ask is whether these 
stars are part of a bound or unbound system, i.e., is 
this an intruding galaxy, or the extended debris of a 
galaxy/star cluster? While information on more stars 
(e.g. giants and horizontal branch stars), would be im- 
portant to address this question fully, some simple esti- 
mates are useful in providing guidance. 

If the system is assumed to be bou nd, we can estimate 
the mass usi ng a virial estimator fe.g.. lHeisler et al.]|1985t 
IGouldl[l993h : 



2G {1/ Rj_^ij)i^j ^ ^ 

where a is the velocity dispersion of the system, R±^ij is 
the projected distance between each two stars. In order 
to measure the velocity dispersion we must account for 
the errors in our velocities. There are two main (identi- 
fiable) sources of velocity error. The primary source of 
uncertainty is the error due to the signal-to-noise ratio, 
which is reported in Table 1. A secondary source of un- 
certainty comes from our imprecise knowledge of the zero 
phase point due to the finite sampling of the Stripe-82 
data. This error is k/{V^N) = 65kms~^/A^ where N is 
the number of data points in each light curve. For our 
data, this has a maximum value of 1 km and therefore 
does not affect our dispersions significantly. The phase 
errors induced by the period errors are similar. 

Assuming all 5 stars in the main velocity structure 
are part of the same physical structure, we obtain a ve- 
locity dispersion of G.Sl^j'g kms^^ which yields a mass 
of i.4tJ:^ X io^Mq. a more conservative possibility is 
that only the "clump" of four stars near (RA,Dec) = 
(356°,— 0.3°) is part of a bound structure, from which 
the star at ^ 352° either has been tidally stripped or is 
not physically associated with the main structure. In this 
case. Equation ^ yields S.Stg^i^ x lO'^M©, whose la 
range is not far from the typical va lue of 1 x 10^ Mp^ for 
the mass interior to 300 pc found bv lStrigari et al.l ()2008D 
for dwarfs with a wide range of luminosities. Moreover, 
the apparent 0.5° radius of the concentration corresponds 
to about 750 pc, rather than 300 pc, and this difference 
could account for the level of discrepancy in mass. Hence, 
this structure is possibly virialized, and if so represents a 
new satellite of the Milky Way. Given the small number 
of confirmed stars and the large physical size, however, 
we cannot at this point rule out the possibility that it is 
being (or has been) tidally disrupted. 

If the apparent structure is physical (and not a chance 
superposition) then whether it is bound or not, the pres- 
ence of 4 RR Lyrae stars within a ~ 0.5 deg^ rectan- 
gle implies that there must be a concentration of other 
stars within the system. By the fuel consumption theo- 
rem, Srr, the RR Lyrae specific frequency (the number 



per y-band lumino sity, normalized to Mvnorm = —7.5; 
ISuntzeff et aI1[T99ll) can be cast in terms of the fraction 
of He burning that takes place within RR Lyrae stars, 
VRR- 
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0.4AMboi 



360r]RR (2) 



4mH - mHe 1 - 7?Ho 
where mx is the atomic mass of X, ryjjo ^ 25% is the ini- 
tial abundance of He, AMboi = Mv,^^ — Afv',norm + ABC, 
Mv,RR = 0.6, and ABC = 0.45 is the difference in bolo- 
metric corrections between RR Lyrae stars and typical 
iant stars. The highest observed value is Srr = 158 
Siegel et al.ll2001[ ). i.e., r]RR = 0.44. The surface bright- 
ness of the Pisces concentration may be expressed as 
V = (34.0 -I- 2.51og77fli^) magarcsec"^, and hence to be 
"recognizable" a.s a V = 30 magarcsec"^ overdensity 
would re quire 7?/?^ < 0.025 or Srr < 9. Such systems 
do exist (jHarris et al.lll996l ). but they are far from uni- 
versal. We cannot therefore, probe whether the Pisces 
structure is a physical (or chance) association based on 
the presence (or absence) of an already known photo- 
metric counterpart in giants. Deeper imaging data of 
this region as well as further spectroscopic follow-up of 
giant and horizontal branch stars in this region of the 
sky is necessary to confront this hypothesis. 

4.3. Unbound Case: Galaxy or Globular Cluster? 

If the object is an unbound stream, as opposed to a 
bound system, we must ask whether this represents the 
disrupted remnant of a globular cluster system or the 
remnant of a satellite galaxy. Robust characterization 
of streams and disrupting galaxies, both theoretically 
and observationally, is a challenge. The predicted phase- 
space morphology depends on a variety of initial con- 
ditions including initial mass, M/L ratio, stellar density 
profile and orbit parameters, and the observations are de- 
manding and not always conducive to multiplexing due 
to the low density of targets. A thin physical extent 
transverse to the direction of motion on the sky, uniform 
stellar color-magnitude diagram and narrow velocity dis- 
persion would typically indicate a globular cluster stream 
as opposed to a galaxy stream. The large width of our 
structure suggests that if it is disrupting, it is a satellite 
galaxy as opposed to a globular cluster. The velocity 
dispersion is comparable to what is observed within the 
Sagittarius stream (Majewski e t al. 2004) or with in the 
Anticenter Stream ((Grillmair. Carlin fc Maiewskil l2008') . 

If the Pisces Overdensity is indeed a disrupted satel- 
lite it is interesting that located near the central core 
of this structure are three stars with similar velocities 
but offset from the main structure by nearly 100 km s^^. 
As such large velocity offsets are not plausible simply by 
disruption effects alone, these stars are either a separate, 
perhaps bound, structure or they are random interlopers. 
The latter scenario is of little interest. In the former sce- 
nario it is intriguing, but perhaps not surprising, that we 
observe two overlapping streams at this position in the 
halo. More locally, the Sagittarius system is so extensive 
and has made a sufficiently large number of orbits that 
many objects have been found within it (e.g. Segue 1) 
that can only be disentangled using velocity and metal- 
licity information. Farther out in the Galaxy, Bell et 
al. (2008) have determined that at least 50% of the stel- 
lar density is in a clumpy form. Detailed comparison 
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with cosmologically-motivated models for the spaghetti- 
like nature of the halo would be useful in determining 
whether such self-overlapping systems are rare or com- 
monplace in the Milky Way at these distances. 

5. CONCLUSIONS 

We present spectroscopic confirmation of the photo- 
metric overdensity observed in RR Lyrae stars toward 
the constellation Pisces. We suggest that this system is 
a dwarf galaxy - possibly in the process of disruption, 
possibly already disrupted, or possibly bound with very 
low surface brightness. 

The nature of the stellar halo at distances of ^ 100 kpc 
is still relatively uncharted territory. Using RR Lyrae 
stars to explore this regime of phase space is a novel and 
exciting way forward. Large photometric datasets such 



as those provided by SDSS and similar experiments in 
conjunction with large ground-based telescopes to obtain 
spectra should yield important insight into the nature of 
the Milky Way's perhaps troubled, perhaps tranquil past. 
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TABLE 1 
RR Lyrae Targets 
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Periods derived in g and r-bands were identical for all but 1 star. In that case the velocity changes from -197.7 kms 
^ . Quoted velocities are heliocentric and values in parenthesis are the estimated observational errors. 
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Fig. 1. — RR Lyrae stars in Stripe-82. Top Panel; positions of RR Lyrae stars with magnitudes 19.9 < Vo < 20.8 in the SDSS Stripe-82, 
color coded in 0.1 mag intervals, brightest to faintest; (gray, red, orange, yellow, green, cyan, blue, magenta, purple). Notice the extreme 
concentration at RA —25° — 0°. Bottom Panel; Zoom of Panel A in the region of the concentration. 
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Fig. 2. — Lightcurves for targeted RR Lyrae stars in Stripe-82. Lightcurves are ordered as they appear in Table 1. The top five are part 
of the common- velocity structure near — 75kms~^. 
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Fig. 3. — Histogram of heliocentric radial velocities for 8 RR Lyrae stars in the Pisces Overdensity. The concentration at — 75kms~^ is 
obvious. The clump at ~ — 190kms~^ is also of note. 
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Fig. 4. — RR Lyrae targets for this program. Green, blue, cyan symbols correspond to median magnitude ranges of 20.5 < gmed < 20.6, 
20.6 < Qmed < 20.7, Qmed > 20.7. Star symbols correspond to targets with velocity greater than -90kms~^ and triangles are all other 
targets for which we have obtained a radial velocity. Squares show objects for which we do not yet have spectra. 



